Simple microscopic model of the magneto-electric effect in non-collinear magnets 
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A 2-site 2-electron model with s- and p-states on each site, having constrained average site-spins 
S a ,Sb, with angle between them, simpler than the previous closely-related model of Katsura et 
al (KNB), is considered. Intra-site Coulomb repulsion and inter-site spin-orbit coupling Vso are 
included. The ground state has an electric dipole moment it consistent with the result of KNB, 
i cc Rx (S a x St) = f , R connecting the two sites, and application to a spiral leads to ferroelectric 
polarization with direction f , also in agreement with the previous result. However, the present result 
shows the expected behavior n — > as Vso — > 0, unlike the previous one. 

PACS numbers: 75.80+q. 71.27.+a. 77.80.-e. 71.70Ej 



5h 

43 



c 

o 
o 



> 

oo 
o 

o 

-I— > 



c 

o 
o 



X 



INTRODUCTION 

The magneto-electric (ME) effect, i.e. the production 
of an electric polarization (magnetic moment) by appli- 
cation of a magnetic (electric) field, has been known for a 
long time. Multiferroics, e.g. materials that show ferro- 
magnetism and ferroelectricity coexisting, have also been 
known where spins and electric dipoles order at different 
temperatures. Quite recently a new class of multifer- 
roics has been found, in which the magnetism and fer- 
roelectricity are strongly coupled; e.g., they order at the 
same temperature. For this effect, however, the mag- 
netic structure cannot be simple ferromagnetism or a 
collinear spin state; the (vector) magnetization density 
must show spatial variation in its direction as, e.g. for a 
spiral spin state. I refer to the introductory paragraphs 
of several recent papers for a more detailed history and 
list of references. Q IS S Until very recently, the 
effect has been found in antiferromagnets, e.g., simple 
spirals, with no net magnetic moment. In 5], however, 
the magnetic structure of the material studied, C0CT2O4, 
is approximately a ferrimagnetic spiral (there is a net mo- 
ment). 0,01113 

The work of Mostovoy is a phenomenological the- 
ory (see also H3)> whereas microscopic models exhibiting 
the effect are presented in Q] and 0]. Both of the lat- 
ter involve superexchange (where electron hopping be- 
tween magnetic ions involves an intervening oxygen ion). 
In Q the situation is considered where a nearest neigh- 
bor pair of magnetic ions has an inversion center at the 
mid-point, so that there is no Dzyaloshinsky-Moriya in- 
teraction (DMI) [H02; 

the polarization comes from a 
distortion of electronic density without ionic or atomic 
displacements; also the t^ g orbitals considered on each 
magnetic site are chosen to diagonalize the spin-orbit 
coupling (intra-atomic SO coupling), partially removing 
degeneracy of these orbitals. In contrast, the essential 
mechanism invoked in Q depends on the existence of the 
DMI plus electron-lattice interactions and orbital degen- 
eracy (Jahn- Teller effect), and the polarization results 
from ionic displacements. 



It seems sensible to ask, is it necessary to consider such 
complicated models to obtain an essential understanding 
of this unusual ME effect? We present here a much sim- 
pler model which yields the effect, namely the electric 
dipole moment for a pair of magnetic sites, Q 
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(1) 



where R connects the sites a and b, and S a ,S(, are the 
average spins at the sites. It is most closely related to 
that of 3], in that it considers the interaction of two 
magnetic sites with a center of symmetry (so there is no 
DMI). Also, the source of the spatial variation of the or- 
dered spin density is due to outside effects of spin-spin ex- 
change interactions, as in 0]. The hopping is direct-there 
is no oxygen, and no ionic displacements or orbital de- 
generacies. It involves inter-atomic spin-orbit coupling, a 
mechanism different from the model of |3J ; it is obviously 
different from that of [jj] . We add that these characteris- 
tics are more appropriate than the others for C0CT2O4, 
where there is no orbital degeneracy. 

A further motivation is a difficulty with the work 
of . In the 3d transition metal ions, the spin-orbit cou- 
pling V so is the smallest of the various energies, namely 
Coulomb interactions, transfer or hopping integrals, and 
(cubic) crystal field splitting. But in [3j, V so is taken 
as essentially infinite: it is diagonalized, along with the 
crystal field, before the hopping energies are considered. 
This leads to a doublet (T?) and a higher energy quartet 
(Tg), and the quartet is dropped. This is probably the 
reason why tt does not show the expected vanishing when 
V so — > 0. Our model does not have this difficulty. 

Our simple model has two atoms or ions, a and b, ly- 
ing on the x-axis; each has an s-type orbital and 3 p-type 
orbitals lying an energy A higher. Also, there are two 
electrons, I per site. Later we describe briefly the gener- 
alization to the case where each site has 3 ti g electrons, as 
for Cr 3+ on an octahedral site, appropriate to the B-sites 
in XCr 2 4 , X=Co,Mn. 

For our simple model we first constrain the 1-electron 
basis set (including spin) to make the average spins < 
S a > , < > lie at particular angles in the x-y plane 
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and make an angle between them. The 1-electron basis 
states are 

s a = u(r-K a )xa, s b = u(r - TL b )xb 

Pa = w(r-n a )xa, pf = vv(r-R a )xb, (2) 

where the spin states are 

Xa = (a + e^°l3)/V2 

Xb = {a + e^p)/V2. (3) 

(a,/3 are the usual "up, down" states along the z- 
direction.) And we take <fib — 4>o + © as in FIG. 1. 
This spin arrangement is chosen for simplicity. Also, 




FIG. 1: General spin configuration in x-y plane. 

v = x, y, z, R a and R;, are the locations of atoms a and 6, 
and the origin is at the mid-point between the two atoms. 
Finally the "radial" functions u, v and the one-electron 
potential V(r) are assumed to be invariant under y — ► — y 
and z — > —z, V also being even in x. We can allow hop- 
ping for any orbital to any orbital, although, again for 
simplicity, we assume no intra-site transition matrix ele- 
ments, and that the intra-orbital Coulomb repulsion is so 
large as to exclude such inter-site hopping. We include 
also the intra-site inter-orbital (s-p) Coulomb repulsion 
Uq. Finally we include the essential hopping processes 
Sa Pb i s b Pa caused by the spin-orbit interaction 

V SO = c VV x p • s, (4) 

with c = H 2 /(2m 2 c 2 ) : (p, s) = (momentum/?!, spin/7l). 
The crucial matrix elements are, e.g., 

< Pa\V SO \sb >= C < UV a \VV X p|u & > • < XaMXb > . 

(5) 

Fve put v a = v(r — R a ), etc. It is convenient to consider 
the spatial factor here, 



=< w a \(VV x p)j\u b >, 



(6) 



where j = x,y,z. Using the symmetry properties of 
u, v, V stated above, one can see readily that all the quan- 
tities m v { vanish except m v z and m z y . The required spin 



factors in JSJ are straightforwardly found to be 

<Xa\s z \ X b> = (l-e 40 )/4^ 2 

< Xa|* w |X6 > = \[e- i4 "> - e l ^ +e )] = i v . (7) 

We will also need 



77 =< Xa\Xb >= -(l + e m ). 



(8) 



Consider first the simplest case, (fio = —0/2, which 
implies £ y = 0. Then we need only the term m|, which 
is seen to be 



ml = 17, 



(9) 



where 



f A , , fdV d dV d\ , , , , 
^ = - C °J drV ^ )y {^Yy~^) Ub{r) (10) 

is real; further, there is no symmetry reason for this to 
vanish. Thus the remaining SO matrix element JSJ) is 



M = 



11) 



The basic symmetry of the situation allows the as- 
sumption v a Vb,u a <-> Ub under x — > —x. It follows 
that m v z — > —m v z under x — > —x, i.e. under a «-> b. 
The above results plus hermiticity of W x p give all the 
relevant matrix elements. 

Because only the p v orbitals are connected to the 
ground state orbitals s a ,Sb, we can drop the other p- 
states. This leaves four 1-electron states, © with v = y, 
and therefore six 2-electron states. We write them con- 
veniently in terms of j4j,^4 p f, which respectively, create 
an electron in states s a ,Pa', similarly for B\, etc.: 

$ a = A\Bi\0>, $ 2 = At J Bt| >, $ 3 = AlBt\0>, 
$ 4 = 4f?t| >, $ 5 = A\Al\0 >, $ 6 = £?|f?t| ^L2) 

Our final simplification before writing down the Hamil- 
tonian H is the Hubbard-like assumption where hop- 
ping is a 1-electron operator, the essential contribution 
from the Coulomb terms being the intra-site inter-orbital 
Coulomb term Uq. Thus, in the basis ($i, • • • $6), 
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(13) 



The real quantities t, t' are ordinary (kinetic + Coulomb) 
hopping terms: t hops s a to s&, t' hops p a to pb- [13 

Let us calculate the electric dipole moment in pertur- 
bation theory, where all hoppings are small. A non-zero 
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value occurs in first order, surprising to us, having ex- 
pected the dipole moment to come from intra-atomic 
mixing of s-like and p-like orbitals (but this occurs in 
higher order only). The ground state wave function to 1st 
order can be picked off from the first column of H, (|13|l . 
and using the unperturbed energy of the doubly-occupied 
sites: 



*9 = *1 



A + C/o 



(14) 



The dipole moment operator is tt = e(rj + r^). To 1st 
order, one needs the results 

< $i|tt|$5 >= - < $i|tt|$ 6 >*= yy < X b\Xa >, (15) 

where 



V = J d 3 r u(r-R a )y 2 v(r-R b ). 



(16) 



Note that y has dimensions of length. The x and z com- 
ponents of tt vanish by symmetry. 

Then, after a bit of arithmetic, we find 



< TT > = < * g |7r|^ g >= -y 



U + A : 



i®. 



(17) 



Our spins being in the x-y plane, this result is clearly 
consistent with 0J. 

It is instructive to consider the electron density, 



/x , \i , x2 7 sin® 
n(r) = u a {r) + u b (r) 



2(U + A) 
where the "overlap density" is 

n ov (r) = u a (r)v b (r) + u b (r)v a (r) 



y n ov (r), (18) 



(19) 



Thus the charge density responsible for the dipole mo- 
ment exists mainly between the two sites. Note that the 
result — > as Vso — * 0, as expected. 

We now generalize to the case where S and S& respec- 
tively make angles 0o and (f>o + ® with the positive x-axis, 
and remain in the x-y plane (see Fig.QJ. We will see that 
the dipole moment is independent of 4>o . 

But the non- vanishing of £ y forces consideration of p z a 
and p b , due to 10 and the above statement that m y ^ 0. 
Thus the number of necessary 1-electron states is in- 
creased from 4 to 6, increasing the number of 2-electron 
states to 15. Nevertheless, it is again quite simple to 
write down the ground state to 1st order in the hopping 
terms, yielding the generalization of l|14|l . which in turn 
gives the electron density: 

n ( r 1 

n(r) = u a {rf+u b (v) 2 +2 T ^^[ 1 Im^ z n*)y+ 1 l Iin^ y n* 



But from the definitions Q and JSJl, 
Im{£ v rf) = 0. 



(20) 
(21) 



Hence the term cx z doesn't enter, and so 7' is irrelevant, 
and the result reduces to the previous expression (|TH|) . 

To finish this consideration of the 2-site model, we con- 
sider the case where the two spins lie in the y-z plane, 
i.e. perpendicular to R. Then the spin states are 

= [cos{9j2)a + tstn(9^/2)/3]/V2, /1 = a, b. (22) 

In this case one can easily see that £ z ,£y,?? are all real, 
so that H2()|l says the dipole moment vanishes, again con- 
sistent with (|T|). 

Now apply these results to a crystal where the spins 
form a spiral [lit ITil Hi| 

< S„ >= l/2[xcos(k • R„) + ycos(k ■ R„)] (23) 

For simplicity we can consider a cubic crystal and the 
propagation vector k along a principle cubic direction, 
say x. Then k • R„ can be written n®, where Q is the 
spiral turn angle. By choosing <fio — n®, one sees from 
FIG.nthat a spiral is generated by increasing n by steps 
of unity. Then the induced change in electron density in 
the bond j, j + 1 is 



Sn(r] 



3,3+1 



A) 



3,3+1 



2{U + A) 



^ysin®. 



(24) 



That is, the dipole moment in each bond is the same 
(ferroelectric ordering) , and the overlap density is merely 
the translation of this density from one bond to the next. 

This way of generating results for a crystal from those 
of a pair of magnetic atoms follows that of Q . 

In summary, the present microscopic model yields an 
electric dipole moment tt resulting from canted spins, in 
the direction of f = R x (S a x S b ) with the same depen- 
dence, sin®, on the angle Q between the spins. When 
applied to a crystal with a simple spiral spin state, the 
f component of tt is the same for each bond, yielding a 
ferroelectric state, just as in the previous closely related 
theory of Katsura et al 0- In contrast, however, our 
result — > as Vso — * 0. 

We mention a very rough estimate of the size of the 
effect. We assume hydrogen Is and 2p orbitals, R=3A 
(the Cr-Cr distance in CoCr 2 4 ), U +A=2eV, and V the 
potential energy of an electron in the field of 2 protons. 
We find the coefficient of sin® in l(T7|l (the dipole moment 
for one bond), to be 2 x 10~ 36 C-m (Coulomb-meters). 
Assuming a simple cubic lattice with spiral propagation 
along a principle cubic direction, and a volume per site 
of (3A) 3 , yields a polarization of ~ 0.1 /iC/m 2 , about an 
order of magnitude smaller than found [3 in CoCi^O^. 
In view of the crudeness of this estimate, it suggests that 
2] -the mechanism is probably relevant to real materials. 

Finally, we extended our model to the case where each 
site is like Cr 3+ in an octahedral field (B-site), the three 
3d-electrons being in the high-spin state i| g j as m Co 
and Mn chromite. The p-states are from the 4p shell. 
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Also note that the dominant B-B exchange interaction is 
direct, the superexchange being 90° |l6j, suggesting that 
our simple model neglecting an intervening oxygen might 
be somewhat realistic for the B-B pairs in these spinels. 
Our motivation is to check that our basic mechanism is 
robust in going to a more realistic model. 

The model Hilbert space is now much larger. Our basis 
functions are constructed as follows. The 3 t 2g orbitals 
are t v (r) = g u (r)u(r),gi = xy,g 2 = yz,g 3 = zx. Sim- 
ilarly the p-orbitals are p v (r) = h u v(v),hi = x, h% — 
y, h 3 = z; u(r) and v(r) are invariant under inversion, 
r — > — r. The full 1-electron basis states are 



plus 3 terms for o x . But in calculating it, these ma- 
trix elements get multiplied by corresponding elements 
of the position or displacement operator as in r^tVy „, 
and many of the position matrix elements vanish by sym- 
metry. It turns out that only the o z terms contribute, and 
the final result for the dipole moment is 



< 7T > = 



Uo 



-sm9(7ij/i +722/2), (31) 



where 



rpU 
b 

b 



tf(r) X . 

b b 

K(r)x ? ; 

b b 



(25) 



£„(r) = t v (r — R a ), etc. We again simplify somewhat by 
choosing the spin states \ a as m © with (j>o = —0/2. 



Let Ay,Ap create T^,P^ respectively, and similarly for 
Bp , Bp . Write these as C\,-y = l,--- ,12, with 7 = 
1, • • • ,6 corresponding to the T-states, 7 = 7, • • • , 12 to 
the P-states. 

Then the spin-orbit interaction is conveniently written 



Vso = c ^2 Vi^xp, i -s. i 



"£<-y\v s °(r,s)\ 7 '>C\C Y , 



(26) 

where v so (r, s) — coS/V x p • s. Then the 6-electron basis 
states (single determinants) are 



$1 - n?ct| >= |o), 

3y 7 = Cl, C 7 |0), 7 <6, 7 '>6. 



(27) 



For the hole vacuum |0), the C\ are ordered C\C\ ■■■C\. 
The ground state to 1st order in the overlap is then 



$y 7 < 7'|w so |7 > 



7<6,7'>6 



(28) 



where the E are the unperturbed energies. The terms 
that contribute to the electric dipole moment are the in- 
teratomic elements, e.g. 



C <K|Wx P |i£>-<Xa|s|Xfe> 



(29) 



Most of the matrix elements of o vanish by symmetry; 
noting that the y-component of the spin matrix element 
t; y vanishes (by the choice of cj) , as before), we need only 
the components o l , i = x, z. The only ones of these ma- 
trix elements that survive are 



io\ v = c <p 1 a \{VV x p),|tj > 
io z 32 , = c <p 3 a \{X7Vx p) z \t 2 b > 



(30) 



7i 
Vi 

1)2 



<pl\o z \t\>, l2 =<pl\o*\tl> 

dhplytl = (1/2) / d 3 r(x 2 - R 2 /4)n a 



d'rpiyti = (1/2) / d'ry'z'n 



(32) 



Thus the result is similar to the previous one (|17fl . The 
terms, ycii-,% — 1,2 do not vanish by symmetry. Thus 
we have confirmed that the basic mechanism, involving 
inter-atomic SO coupling, gives the interesting ME effect 
in this rather realistic model. 
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